Structural response to local charge order in underdoped but superconducting 

La2-x(Sr,Ba)xCu04 
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We report an anomalous local structural response in the Cu02 planes associated with the 
appearance of charge inhomogeneities at low temperature in underdoped but superconducting 
La2-a:(Sr,Ba)a;Cu04. We used pair distribution function analysis of neutron powder diffraction 
data. The increase in the Cu-O bond length distribution at low temperature has an onset temper- 
ature which correlates with observations of charge and spin freezing seen by other probes. 
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Two unusual phenomena are observed in the normal 
state in the underdoped cuprates: a pseudo-gap [y in 
the electronic [||-g and magnetic |^,0| densities of states, 
and the possibility that the charge density in the super- 
conducting planes of these materials is inhomogencously 
distributed, possibly in a striped morphology ||q-[l8||. It 
is important to establish the role that charge inhomo- 
geneities have in the high-Tc phenomenon itself. Un- 
like the pseudogap phenomenon, their universal observa- 
tion among different high-Tc systems has not been estab- 
lished. The strongest evidence for them in the cuprates 
is the observation of long range ordered static charge 
stripes in La2-a:-i/NdySra;Cu04 compounds |^. These 
have been seen in both insulating and superconducting 
compounds but they appear to compete with the super- 
conductivity p9| . On the other hand fluctuating short 
range ordered charge stripes may play an active role in the 
high-Tc phenomenon ]13-0]. They also give a natural 
explanation for the observation of incommensurate spin 
fluctuations which have been seen in La2-a;Sr2;Cu04 B] 
and YBa2Cu306+5 pOf as well as being able to explain 
various other experimental observations |lC|Jl9| -pi^ . It 
is important to establish both the ubiquity of charge 
inhomogeneities in underdoped cuprates and their rela- 
tionship to superconductivity. Here we present diffrac- 
tion evidence that establishes the presence of tempera- 
ture dependent atomic scale structural inhomogeneities 
at low temperature in underdoped but superconduct- 
ing La2-a;(Sr,Ba)a;Cu04 samples. This observation is 
naturally explained by the appearance of charge inho- 
mogeneities at low temperature in these samples. The 
inhomogeneities appear at a temperature which corre- 
lates with spin and charge freezing inferred from trans- 
port [|9), NQR |2l|, and XANES || measurements. 

We used the atomic pair distribution function (PDF) 
analysis |gj] of neutron powder diffraction data to study 
the local structure of La2-a;(Sr,Ba)xCu04. Structural 
distortions coming from charge inhomogeneities appear 
in the PDF as an anomalous broadening of the near- 
est neighbor in- plane Cu-O bond length distribution [G5[ . 



The average in-plane Cu-O bond length shortens on hole 
doping. This is observed experimentally [p6| and is ex- 
pected on the grounds that holes are being doped into 
a (7* antibonding band [p7[ thus stabilizing the bond. 
Charge inhomogeneities imply a coexistence of heavily 
and lightly doped regions of the Cu02 plane. The lattice 
will respond if the charge inhomogeneities are fluctuat- 
ing on phonon time-scales or slower. This will result in 
a distribution of lengths for the in-plane Cu-O bond and 
correspondingly to a broadening of the atomic pair dis- 
tribution. This can be measured directly using the PDF 
analysis of neutron powder diffraction data. The PDF 
technique, which is common in the study of glasses p4[ , 
is equally well applied to crystalline systems where it re- 
veals precise information about the local atomic structure 
going beyond the approximation of crystallinity pq ] . 

Powdered samples of La2-a;(Sr,Ba)a;Cu04 [x = 
0.125,0.15) of ^ 10 g were synthesized using standard 
solid state techniques [p9^,p0|. The samples were char- 
acterized using x-ray diffraction and susceptibility mea- 
surements. The oxygen content was verified by measur- 
ing the c-axis parameter that was found to fall on the 
expected curve for stoichiometric samples P6| . Neutron 
powder diffraction measurements were carried out on the 
High Intensity Powder Diffractometer at the Manuel Lu- 
jan Neutron Scattering Center (MLNSC) at Los Alamos 
National Laboratory and on the Glasses, Liquids and 
Amorphous Diffractometer at the Intense Pulsed Neu- 
tron Source (IPNS) at Argonne National Laboratory. 
The samples were sealed in vanadium tubes with He ex- 
change gas. Data were collected as a function of tem- 
perature from room temperature down to 10 K using a 
closed cycle He refrigerator. Standard corrections |3l| ] 
were made to the raw data, to account for experimen- 
tal effects such as sample absorption and multiple scat- 
tering, using the program PDFgetN [|2[, to obtain the 
total scattering structure function, S{Q). This contains 
both Bragg and diffuse scattering and therefore informa- 
tion about atomic correlations on all length scales. The 
PDF, G{r), is obtained by a Fourier transformation ac- 




FIG. 1. (a) Reduced total scattering structure function, 
Q[S{Q) - 1], from Lai.g75Sro.i25Cu04 at 300 K. (b) The re- 
sulting PDF, G{r) (open circles). The solid line is a fit to the 
data of the crystallographic model with the difference curve 
below. Arrows indicate the PDF peaks at r = 1.9 A and 
r = 7.2 A whose widths are plotted in Fig. bl 

cording to G{r) = ^ /g°° Q[S{Q) - 1] sinQr dQ, where Q 
is the magnitude of the scattering vector. The PDF gives 
the probabihty of finding an atom at a distance r away 
from another atom. The PDF from Lai.875Sro.i25Cu04 
at 300 K is shown in Fig. H^ (b) with the diffraction data 
in the form of Q[S{Q) — 1] in Fig. ^ (a). Superimposed 
on the PDF is a fit to the data of the average struc- 
ture model using the profile fitting least-squares regres- 
sion program, PDFFIT |||. The S{Q) data were termi- 
nated at Qmax — 23 A^^. This is a conservative value 
for Qmax in typical PDF measurements. The data from 
high-Q has a poorer signal-to-noise ratio because of the 
effect of the Debye- Waller factor. By eliminating it from 
the Fourier transform we improve the signal-to-noise ra- 
tio of our data and the temperature to temperature re- 
producibihty of the PDFs. This reduces the possibility 
that observed effects are noise artifacts. We can there- 
fore have confidence that any effects that survive this 
conservative approach to Fourier transforming the data 
are real. 

In La2_a;(Sr,Ba)a;Cu04, the first peak in the PDF at 
r = 1.9 A originates from the in-plane Cu-0 bond. The 
width of this peak comes from the relative motion of 
nearest neighbor in-plane Cu-0 pairs, plus any static or 
quasistatic bond-length distribution, averaged over the 
whole sample. We have studied the mean square width, 
a^ oc {u^), of this peak as a function of temperature for 
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FIG. 2. Mean-square width, a , of the in-plane Cu-O 
PDF peak at 1.9 A as a function of temperature for 
(a) Lai.g75Sro.i25Cu04, (b) Lai.85Bao.i5Cu04, and (c) 
Lai.85Sro.i5Cu04. The solid line gives the T-dependence pre- 
dicted by the Einstein model [M. The insets show the in- 
verse squared peak height for the same peak with an Ein- 
stein curve superimposed, (d) Temperature dependence of 
l/h^ for the PDF peak at 7.2 A for Lai.875Sro.i25Cu04 (tri- 
angles), Lai.85Bao.i5Cu04 (diamonds), and Lai.g5Sro.i5Cu04 
(squares) Inset shows the strontium-dopin g d ependence of 
a'^{x) for La2-a;Sra;Cu04 at 10 K from Ref. [E5|. 



a series of underdoped La2-x(Sr,Ba)2;Cu04 compounds. 
The results are reproduced in Fig. g. Peak profiles in 
the PDF are well modelled using a Gaussian function 
convoluted with a Sine function, sin Qmaxr/QmaxT [ p5[ . 
The Sine function results from Fourier transforming the 
finite-range data. Since Qmax is a known experimental 
parameter it is possible to extract intrinsic peak widths 
for the underlying Gaussian distributions. The results of 
this convoluted fitting process arc shown in Fig. 0(a)-(c). 
It is clear from the Figure that the peak width decreases 
with decreasing temperature as expected. However, be- 
low a certain temperature the Cu-0 bond length distribu- 
tion broadens on further decrease of temperature. This 
effect cannot be explained by normal thermal behavior 
as indicated by the solid lines which have the expected 
Einstein form |3J]. There also are no structural phase 
transitions occuring at these temperatures. 

The same qualitative result was obtained from the data 
directly without carrying out a convoluted fit. First, we 
simply plot the inverse-squared PDF peak height, l//i^, 
obtained directly from the data. This is a model inde- 
pendent measure of a^ since the integrated area under 



the PDF peaks is conserved |3^. The inverse-squared 
peak heights are shown in the insets to Fig. H(a)-(c). We 
also fitted the 1.9 A PDF peak with pure Gaussian func- 
tions that were not convoluted with Sine functions (not 
shown). Both these approaches reproduced the quali- 
tative result shown in Fig. ^(a)-(c) giving us confidence 
that it has a real origin and is not an artifact of the con- 
voluted fitting procedure. All of these measures of the 
PDF peak width confirm the observation in the convo- 
luted peak fits that the underlying in-plane Cu-0 pair 
distribution is getting broader with decreasing tempera- 
ture below some temperature, Tgi- 

Peaks not involving in-plane Cu-0 pairs, at higher val- 
ues of r, in the PDF behave canonically. This is shown 
in Fig. |(d) where l/h'^{T) of the peak at r = 7.2 A 
(indicated with an arrow in Fig. yj) from each of the sam- 
ples is plotted with an Einstein curve superimposed. As 
expected, no upturn is observed at low temperature. 

The broadening of the r = 1.9 A PDF peak at low 
temperature can be explained if charge inhomogeneities, 
such as charge stripes, are manifesting themselves in the 
structure at low temperature. This will occur both if the 
electronic correlations are appearing at low temperature 
or if preexisting correlations are slowing down and be- 
ginning to interact with the lattice. It was shown in an 
earlier PDF study [g5[ that a gradual broadening with 
increasing doping at 10 K of the r — 1.9 A PDF peak 
in La2-xSr2;Cu04 could be well explained as a micro- 
scopic coexistence of heavily doped and undoped regions 
of the copper-oxygen plane. The 2:-dependence of this 
PDF peak width measured at 10 K is reproduced in the 
inset to Fig. 0(d). This can be compared with the intrin- 
sic peak widths at low temperature from this study. 

The original x-dependent data were interpreted as fol- 
lows. The relatively sharp peaks in the a; = 0, 0.25, 
and 0.30 data were assumed to have a single valued bond 
length broadened by thermal and zero point motion re- 
sulting in a mean-square width of ^ 30 pm^. The rel- 
atively broader peaks observed in the underdoped com- 
pounds {x = 0.05, 0.10, 0.125, 0.15) could be explained 
as a superposition of sharp peaks that are shifted in po- 
sition originating, respectively, from less doped and more 
heavily doped regions of the copper oxide plane p5| . This 
very simple model independent analysis is likely to be an 
oversimplification of the real situation where local strains 
may lead to broader distributions of the PDF peaks; how- 
ever, it establishes unequivocably that the observed ef- 
fects in the PDF are consistent with structural distor- 
tions originating from charge inhomogeneities. Despite 
the current measurements being made on different ma- 
terials at different times using different diffractometers 
it is clear that both the low-temperature thermal width 
of 25-31 pm'^ extrapolated from the Einstein model, as 
well as the excess peak height of ^ 10 — 15 pm-^, are in 
rather good agreement with our earlier observation of the 
x-dependence of La2-a;Sr2;Cu04 p5[. This indicates that 



the underlying origin of the peak broadening is the same. 

The in-plane Cu-0 pair correlation has been stud- 
ied in a number of XAFS measurements [p6[-p8[. The 
data of Lanzara et al. p7[ qualitatively suggest an up- 
turn in the width of the distribution at low temperature. 
However, this result may not be significant since later 
work suggests that uncertainties in unpolarized XAFS 
measurements are larger than the observed effects pS] 
and that polarized XAFS measurements are necessary 
to obtain higher precision pO]. In particular, this latter 
study puts an upper limit of 0.017 A on possible non- 
thermal disorder amplitude present in the in-plane Cu-0 
bond distribution of Lai.875Bao.i25Cu04. This is not far 
from our suggestion of a ~ 0.02 A splitting observed in 
La2-2;Sr3;Cu04 [ p5[ and in the current work. Our data 
will be compared with the result of Haskel et al. in more 
detail elsewere |39|| . 

We have extracted a temperature, T^i, where the 
structural inhomogeneities set in by taking the differ- 
ence, Acr^, between the observed width and the Ein- 
stein curves plotted in Fig. g. The resulting values for 
Tsi are 125 K for Lai.875Sro.i25Cu04 and 60 K and 
100 K for Lai.85Sro.i5Cu04 and Lai.85Bao.i5Cu04, re- 
spectively. These are shown in Fig. || as solid hexagons. 
The estimated error bars are rather large since the exact 
value of Tsi depends on parameters used in the Einstein 
fits; also our data-sets are somewhat sparse. However, 
they give a temperature scale where the effects of charge 
inhomogeneities first appear in the local structure. 

In Figure H we show a phase diagram for 
La2_a;Sra;Cu04 with T^,; plotted along with Tc and T* 
obtained from the literature pi. Superimposed on this 
diagram are T^/, the onset temperature for NQR signal 
"wipe-out" pl| , T„, the temperature where the devia- 
tion of the normalized resistivity, p/aT, reaches a critical 
value ||l9|l , and Tj; , the temperature where an anomaly is 
seen in XANES data |g^ . All of these characteristic tem- 
peratures have been associated with charge or spin freez- 
ing. It is clear that the Ts^'s obtained from the present 
data correlate quite well with the other measures of spin 
and charge freezing. 

Our results clearly show that the charge inho- 
mogeneities are strongly coupled to the lattice in 
La2-2;(Sr,Ba)2;Cu04 compounds and become pinned by 
the lattice at low temperature. In the absence of Nd 
the pinning is not complete and the charges do not or- 
der over long range, even at x=0.125 in the Sr doped 
compound. Nonetheless, their strong coupling to the lat- 
tice will make them quite immobile. Our measurements 
yield the instantaneous structure and we cannot distin- 
guish whether the inhomogeneities are static or fluctuat- 
ing on phonon time scales or slower. It will be interest- 
ing to see whether similar effects are seen in the PDF of 
HgBa2Cu04+5 which is a single layer cuprate supercon- 
ductor like the La2-2;(Sr,Ba)3;Cu04 compounds but has 
a much higher Tc. It is possible that electronically driven 




0.05 0.1 0.15 0.2 

strontium content, x 

FIG. 3. Phase diagram of La2-i:Sri:Cu04 showing the tem- 
peratures of pseudogap openin g , T* |l| , XANES anomaly, 
T^ H, NQR spin freezing, T^/ [|j transport upturn, T„ jl| 
and the Tsi's determined from the present measurements. 
Tsi is known to be below 10 K for La2-a:Sri:Cu04 with 
X > 0.2 psl as indicated. Tc, is shown as solid circles joined by 
a line. The inset is the same phase diagram on an expanded 
temperature scale. 



stripes are important for superconductivity but a strong 
coupling to the lattice is destructive. However, phonon 
anomalies have been associated with charge stripe forma- 
tion in YBa2Cu306+5 [|l^,^ and theories exist in which 
the charge stripes are stabilized by the lattice. Resolving 
this issue will be a key component in gaining a complete 
understanding of high temperature superconductivity. 
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